Reexamination of decoherence in quantum walks on the hypercube 
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The effect of decolierence on tfie continuous-time quantum walk on the hypercube is revisited. 
Previously, an exact solution was found for a decoherence model that preserved the effective tensor- 
product form of the dynamics. Here a new model is presented, inspired by perfect quantum state 
transfer in qubit networks. A perturbative solution is found for the dynamics of this model which 
is not of a tensor-product form. In contrast to previous results, the hitting probability has a lower 
bound that is independent of the hypercube dimension. 
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I. INTRODUCTION 

The quantum walk on the hypercube is a quintessential 
problem in quantum information processing. The main 
uses of quantum walks can be divided into three parts: 
search, hitting, and fast sampling algorithms. The first, 
the search of a database, is based on Grover's algorithm, 
and can be implemented using a quantum walk on a hy- 
percube graph flj . Hitting algorithms are those in which 
a particle can propagate from one node to another node in 
a graph (sometimes called the graph traversal problem) . 
For the hypercube, it has been shown that the quantum 
walk performs this task exponentially faster than clas- 
sical random walks 0; a graph for which the quantum 
walk is exponentially faster than all classical algorithms 
was found by Childs et al. [3]. Finally, a fast sampling 
algorithm approximates a probability distribution (often 
uniform) over some graph starting from a simple initial 
condition — this is quantified by what is called the mix- 
ing time. The hitting and instantaneous mixing times 
for the hypercube were found in the seminal work by 
Moore and Russell [3]. Other important works on the 
hypercube include recent results on the hitting fs'] and 
mixing times jSj jthere have also been general reviews on 
quantum walks 0, [1] ■ 

In addition to potential algorithms, the quantum walk 
on the hypercube can be used to devise qubit networks 
that perfectly transfer a quantum state between nodes 
P|, and could possibly be implemented using supercon- 
ducting qubits 10]. For these physical implementations, 
one must consider the inevitable coupling of each qubit 
to its environment. The simplest such model leads to 
a particular decoherence model, the solution of which is 
the subject of this paper. More generally, decoherence 
in quantum walks may be useful [1, [ll| , and could even 
have implications for biological systems p^ . 

Previous work on decoherence in quantum walks on the 
d-dimensional hypercube began with the work by Kendon 
and Tregenna [llj, who looked at a discrete-time quan- 
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tum walk interrupted by measurements. They found nu- 
merically that by measuring the position of the walk with 
probability p after each time step the first hitting prob- 
ability was exponentially (in d) suppressed. An exact 
solution for decoherence in a continuous-time walk was 
found by Alagic and Russell ^T§\, using a model where 
dephasing occurred between different subspaces of the 
hypercube. Here too, the first hitting probability decays 
exponentially in d. 

In this paper I consider a continuous-time decoherence 
model — here called the vertex model — that is similar to 
the discrete-time model mentioned above [Tlj. The spe- 
cific form of this model arose from a study of decoherence 
in quantum state transfer [13] ■ By analyzing the master 
equation of the vertex model perturbatively, an analyt- 
ical solution is found, good for weak decoherence. For 
this model and a natural choice of parameters the quan- 
tum probability to traverse a d-dimensional hypercube 
remains bounded for all d, in opposition to the exponen- 
tial suppression found previously. 

This paper is organized as follows. The decoherence 
models are presented in Section II along with a review 
of the exact results for the subspace model [1^. In Sec- 
tion HI, a perturbative approach is presented, giving an 
analytical expression for the hitting probability as a func- 
tion of time. In Section IV, I compare this model with 
the discrete and continuous models, and discuss the role 
of information lost due to decoherence. Finally, I con- 
clude in Section V by summarizing these results and their 
implications. An Appendix extends the results of Ref. 
(l(j ] and shows how different models of decoherence in a 
hypercube network of qubits generates the decoherence 
models considered here. 



II. DECOHERENCE MODELS 

In this section I first describe two models of decoher- 
ence: the subspace model studied by Alagic and Russell 
[T3 ] and the vertex model presented here. The quantum 
walk on the hypercube involves the continuous-time dy- 
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namics (with h = 1) 
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where the matrix elements of H are proportional to the 
adjacency matrix of the d-dimensional hypercube. The 
vertices of the hypercube are labeled by bit-strings of 
length d, and two vertices are connected if their labels 
differ by one bit. This induces a tensor-product struc- 
ture, in terms of which the Hamiltonian has the simple 
form 



iJ = 1 ^ ■ ■ ■ cr-r ■ 



1, 



(2) 



where the Pauli matrix occurs at location j. Here and 
in the following, lu and t are taken to be dimensionless 
quantities. Using this Hamiltonian, one has 



l*(0)). 



(3) 



Letting \a) = |0---0) and \b) = the transfer 

probability from |^(0)) = \a) to |6) is given by 

P{t) = \{b\e-'^^a)\^ = I sin(cjt)|2''. (4) 



The first hitting time occurs at T = ■k/{2lo), with unit 
probability. 

In the presence of decoherence, however, this simple 
dynamics is modified. First, the system is no longer de- 
scribed by a wavefunction but rather by a density ma- 
trix. Second, the evolution is no longer generated by 
a Hamiltonian operator, but rather by a superoperator, 
here assumed to be of the general Lindblad form [l^ 



dtp{t) 



i[H, p{t)]+J2 A, {L,pL] - h]L,p- \pL]L, 



(5) 

The Lj are called the Lindblad operators, with rates Aj. 

Alagic and Russell considered Lindblad operators of 
the form: 



1 : 



■H. 
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(6) 



where a is or 1 (these should be included in the sum) 
and Hq = |0)(0| and Hi = are the projectors for co- 



ordinate j. The decoherence rates are all identical \j = A 
(see Section IV for a discussion of these rates) . Note that 
these operators do not project the state onto any partic- 
ular node of the hypercube network, but rather onto one 
of two subspaces in which the node has a (or 1) in the 
j-th bit. For this reason, I will call this the subspace 
model. 

The resulting master equation can be given a tensor- 
product form, which allows for the solution of the density 
matrix by solving for the eigenvalues and eigenvectors of 
a 4 X 4 matrix [13]. Using these and the initial condition 
p{i)) = \a){a\, the hitting probability is 

P,(t) - {b\p{t)\b) 



e 



-At/2 



{cos(/3V2) + |sin(/?t/2)} 



(7) 



with P = Vl6w2 _ x^. 

The decoherence model considered here uses the same 
general form, again with identical rates Xj — A but with 
Lindblad operators of the form 



in 



Hn 



H„ 



(8) 



each a being or 1. These operators project a state 
onto the vertex given by \ai ■ ■ ■ ad)- I call this the vertex 
model, and note that it can be derived from the master 
equation for qubits undergoing dephasing independently 
in a hypercube network, as shown in the Appendix. Note 
also that this decoherence model has essentially nothing 
in common with the graph on which the quantum walk 
occurs. This is in constrast to the subspace model, in 
which decoherence somehow knows the structure of the 
network. The vertex model is somewhat more physical, 
in the sense that for a physical network, decoherence is 
more likely to affect each vertex independently, and less 
likely to know about the full hypercube structure. This 
will be discussed in Sec. IV. 

It is instructive to reexpress the master equations in 
terms of the matrix elements px^y = {x\p\y), where x = 
xi---Xd and Xj = {0,1}. For the vertex model, px,y 
satisfies 



Px,zHz,y) — A(l — Sx,y)Px,; 



(9) 

where Sx.y is the Kronecker delta and Hx^y — {x\H\y). 
This master equation involves the decay of all off- 
diagonal matrix elements with rate A, consistent with 
dephasing of every vertex. The subspace model has 



dtPx,y = X! i^x,^P^,v 



Px,zHz,y)-X} ^{'^-5x,,yj)Px,y, 
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where 6xj,yj is the Kronecker delta for the j-th bit of x 
and y. This too involves the decay of the off-diagonal 
elements, but these are now weighted by the number of 
bits in which the indices x and y differ. Those elements 
with indices differing in n bits decay with a rate nA. This 
is somewhat surprising — the subspace model presumably 
gets less information about the location of the "walker" 
in the quantum walk, yet in the end causes more deco- 
herence than the vertex model. This will be discussed in 
Section IV. 



III. PERTURBATIVE SOLUTION 

In the limit of weak decoherence A ^ cj, perturba- 
tion theory can be used, starting from the eigenstates 
of H . However, since this is a master equation, careful 
treatment of degeneracies is required. These degeneracies 
could be analyzed using an angular momentum repre- 
sentation or perhaps other group theory methods. Here, 
however, I will take a direct approach, requiring the fol- 
lowing notation. The density matrix of the hypercube 
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can be rewritten as a vector, using the basis states 

|xi • • • xd; yi - ■ - yd) ^\xi - ■ ■ xd) (2/1 ■• • yd\, (11) 

where each Xj and yj are either or 1. These states wih 
also be written using an abbreviated notation \x]y). In 
this basis, the superhamiltonian Ti acts as 

'H\xi ■ ■ - Xdiyi - ■ - yd) ^ -i[H, |xi • • • xj) • • • yd\] (12) 

with the value 

T-(-\x;y) = -iYf.^i\xi---xyXd\yi---yd) 



E,=i l^^i ■■■xd]yi---yj- --yd) 



(13) 



where a; = 1 — a; is the bit-flip of x. The Lindblad super- 
operator £0 is defined similarly, and is diagonal in this 
basis: 



^o|a;;y) = -A(l - 5^^y)\x]y). 



(14) 



In addition to this basis, I also introduce the eigen- 
states of 7i, which have the form 

|a;;y).=2-'^^(-l)-'+^-^'|x';2/'), (15) 

x' ,y' 

using the abbreviated notation for the states, with 
x, y, x' , y' representing bit strings of size d, and thus there 
are 2^'^ terms in the sum. Note that I have also used the 
bit-wise product 



X ' X — ^ X J X J . 



(16) 



These states satisfy 



n\x]y)^ ^ 2iuj ^Yy^j ~ yAj \^'^y)x- (17) 

While H is diagonal in this basis, Cq is not. To find 
how it acts in this basis, I first use Eq. and the 

definition of \x;y)x to find 

Co\x; y), = -\\x- y), + A2-'^ ^(-l)--'+^-'|x'; x'). 

x' 

(18) 



Then, using the inverse of Eq. p5)) 

|x;y)=2-'^^(-l)-'+^-^'|x';y')x (19) 

and the fact that 

x' 

where x © x" is bit-wise addition, I find the result 



LQ\x;y)^ = -\\x\y)^^'Xl ^ 4ex" 

x" ,y" 



,y®y"\X ,y / 



Now, since this final sum is over all bit-strings, I can 
make the substitutions x" = x® z and y" = y (B z' and 
sum over z and z' . These substitutions reduce the delta 
fmiction to 8z,z'^ with the final result 

-Cola;; y)x = -\\x; y)x + ^ |x ® z; y ® z)^. (22) 



If all of the states of Ti, were nondegenerate, I could 
simply assert that the effect of decoherence on each 
state is given by {x]y\C{)\x]y)x- However, there are 
many degeneracies, and so degenerate perturbation the- 
ory is required. From Eq. ([TT]) I find that there are 
2(i+ 1 subspaces, with eigenvalues e„ — 2iuj{d~ n), with 
n — 0, • • • , 2d. The degeneracy of subspace n can be 
found to be {2d)\/[n\{2d—n)\\. In the following, I will con- 
sider £0 projected onto these subspaces (with the same 
symbol) and analyze how it mixes those states within 
each subspace. 

Finding the correct form of the eigenstates in these 
degenerate subspaces requires some additional notation. 
Consider, for example, the subspace with n — 2. In this 
subspace, there can be no more than 2 bits of y that are 
unity, and no more than 2 bits of x that are zero. Letting 
those bits occur at positions i and j (with i < j), this 
can be achieved in one of four ways: 



N;j)2 
|j;«>2 



loa- 
10a • 

10,1 



10,1 



1---1:0-- 



•l;0---0). 
•l;0---01jO---0)^ 
••l;0---01,0---0), 
■01,0---01jO---0)^. 



(23) 



Note that the indices i and j label the zeros for x and 
the ones for y; these can be interpreted as the locations 
of "excitations." Using Eq. ((22)) . when projected to the 
n = 2 subspace the Lindblad operator has the following 

effect 

/:o|«)2 = -A|z;)2 + A2-''(|zj; 0)2 |z; i)2 + |j; i)2 + 10; u)2) 

(24) 

where \v)2 is any of the states in Eq. ((23|) . These are for 
i < j ■ For i = j, however, 

/:ob;j)2 = -A|j;j)2 + A2"''^|fc;A:)2. (25) 

k 

Thus, for the states in this subspace the effect of Cq is to 
either mix excitations between the x and y indices, or to 
distribute "paired" excitations over all possible locations. 

States with higher n can be defined and the effect of 
Cq can be analyzed in an analogous fashion. For exam- 
ple \ji • • ■ jn\^)n denotes a state \x]y) with n zeros in x, 
located at positions ^1,^2, • • ■ ,jn, and no ones in y. To 
represent the effect of £0, it is convenient to define the 
symbol 5 as a type of shift operator, transferring the last 
unpaired excitation in x to y, such that 



(21) 



S\jl - ■ ■ jn;^)n = - ■ ■ jn~l]in)n 

^^lil • • •in;0)n = b'l • • • jn~2;i«-ljr, 

S^\n---]n\%)n = |0;ji---j„)„. 



(26) 
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In addition, I define Vm.n to represent a permutation of wfiere 1 = (1, • • • ,1). In fact, combining this result with 
jn with jVn- For example, Eq. ([28|) . and using the identity 



^m,n\il ' ' ' jm ' ' ' — 2 t jn—ljn) n — \jl ' ' ' Jn — 27n; 7mjn— l)n; 

(27) 

where I have reordered the labels to increasing order, 
with m < n — 1. Using these two symbols, I now argue 
that a correct set of states that diagonalize Cq (when 
projected to the relevant subspace) is given by 



|J 



1 ■ • - jn;! 



(28) 

where j = {ji, ■ ■ ■ , j„}, s = (si, • • • , s„) and Sm = ±1- 
When this product is expanded, it produces a superposi- 
tion of states with the indices j distributed in all possible 
ways (2") over the x and y partitions, and with corre- 
sponding weights given by the products of Sm- Note that 
the effect of Cq can be expressed quite simply using this 
basis. Specifically, for \v)n any state formed by shifting 
or permuting |ji • • • j„; 0)„: 



/:ok)n = -A|tO„ + A2-''+"/2|0.j)„, 



(29) 



n(l + ^^™) = 2"%. 



(30) 



the following result is obtained: 



/:o|</',,s-)„ = -A(l-2"-%j)|0,-,-) 



(31) 



This enumerates the states with ji < j2 < • ■ • < jn- 
For cases such as j„_i — jn, where two labels are paired 
(one with an excitation at site j„ for x, the other for y), 
the states have a slightly different form. With p pairs in 
subspace n, the pairs can be distributed in dp = {d ~- n + 
2p)l/[p\{d — n + p)]] different ways. Here I choose to use 
a discrete Fourier transform to construct states for these 
pairs, with the form 



/ n — 2p 



.V = ^n,p n [1 + S,nS]V,n,n^2p ^ e'^.^/C^)/'^- \j, ■ ■ ■ j^_^pk, • • • fcp ; fci • • • fcp) „ , (32) 
\m=l / {k} 



r 



where here S and V act only on the unpaired indices and 
the sum is over all integers with 1 < fci < fc2 < ■ • ■ < 
kp < d and each ^ {ji, • • • ,jn-2p}- The function 
f{k) labels each configuration oi k — {ki, • • • , fc„} from 
to dp — 1 (its explicit form is not needed here), and 
< q < dp — I. The normalization factor is 

«r _ o-n/2+P,-l/2 _ ^-n/2+p ( p\{d - n + p)\ \ ^'^ 

^^'P^^ > \{d-n + 2p)\) ■ 

(33) 

Note that for p pairs there are d\/[{n — 2p)\{d — 7i -I- 2p)!] 
ways of choosing the index set j, 2""^'' possible choices 
of s, and dp possible values of q. Thus, the total number 
of these states is (for n < d) 



ln/2\ 

E 



d\ 



{n-2p)\{d-n + 2p)\ " 



dr,2''-^P = 



(2d)! 



nl{2d~n)V 

(34) 

which — as required — equals the total degeneracy of sub- 
space n. Using an argument similar to the case of p = 
given above, these states can be shown to satisfy 

•^o|0j,s,g)n,p = — A ^1 — 2" ^^dpSg-^Sqfi^ \<t'j.,s,q)n,p- 

(35) 



This completes the analysis of the effect of Cq on each 
subspace. To calculate the hitting probability starting 
from the initial state \a){a\ and ending at I first 

expand these in terms of the eigenstates of Cq on each 
{n,p) subspace: 



■J,s,q — n,p\Yj,s,q 



|a;a) =2"/2-P-'^4/25,.jJ,,o (36) 



and 



Denoting the eigenvalues of |0j,r,g)n,p by 2iu!{d—n) — Xpr, 
with 



\ ^ \ (^ - on-d-2p {d-n + 2py. 

pl{d-n+py. 

the net result for the probability is 



(38) 



_ Y^2d Y^ln/2\ ^.^ „2ioj{d-n)t„~\p„t 

— Z^n=0 l^p=pmin ^j,s.q "'J,s,q'^3,s,qf^ e , 

(39) 

where Pmin = max(0,n — d). Substituting for Cj^g^q and 
dj^s^q and performing the sum over j, s, and q leads to 
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the result 



2d L"/2J ,ir)n-2p-2d( i\d-n 



n— p—Pii 



pi {n — 2p) ! (d — n + ! 



(40) 

This can be simphfied to some degree, by noting there 
is a symmetry between terms with n < d and n > d 
(this can be found expHcitly by letting n' — 2d — n and 
p' ^ p + d — n for n > d). Doing leads to the the final 
result 



Pv{t) = E(-l)''""5nW cos {2Lotid - n)) 



where I have defined the functions 



p=0 



p\{n - 2p)\{d - n + py. 



(41) 



(42) 



IV. COMPARISON OF MODELS 

The subspace and vertex models initially seem some- 
what similar. Each can be written in terms of functions 
that recall damped harmonic oscillators, and indeed, they 
agree for A ^ w and for d = 1: 

P,{t,d=l)Ki Py{t,d=l) = - (^1- e"^*/2 cos{2ujt) 

^ . . . 

For larger values of d, however, significant differences 

emerge. 

For illustration, I compare the two models with uj = I, 
A = 1/5, and for d equal to 1, 4, and 10. In Fig. [1] the 
hitting probability Py{t) is shown for the vertex model, 
using the perturbative solution presented in the previous 
section. Note that this value of A was chosen so that 
the probabilities decayed sufficiently at the transfer time 
T = 7r/(2w) = 7r/2; here AT = tt/IO « 0.314. Direct 
numerical simulations, using a split-operator algorithm 
[m to evolve the density matrix by Eq. ^ show that for 
these values of A and d, the perturbative and numerical 
solutions are in good agreement. The hitting probability 
for the subspace model Ps{t) is shown in Fig. [21 There are 
several common features found in these two figures. First, 
there is clearly damped oscillations of the probability, 
with a frequency approximately given by 2lo. Second, for 
long times the probability becomes constant and equal 
to 2-'^. 

There are, however, some important differences be- 
tween the two models. First, it is clear that the prob- 
ability oscillations decay more quickly in the subspace 
model with d > 1 than in the vertex model. Second, 
while the maximal hitting probability (which occurs near 
t = T = tt/2u} = tt/2) decreases with increasing d for 
both models, it does so much more dramatically for the 
subspace model. This is further explored in Fig. [31 As d 
increases, the probability Ps{T) for the subspace model 



l.Or 

0.8 
0.6 
0.4 
0.2 



ii 



5 10 15 20 

t 



FIG. 1: The hitting probability Pv{t) as a function of time 
for the quantum walk with decoherence in the vertex model. 
The parameters are a; = 1, A = 1/5, and d = 1 (solid), d — 4 
(dashed), and d = 10 (dotted). 
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FIG. 2: The hitting probability Ps{t) as a function of time for 
the quantum walk with decoherence in the subspace model. 
The parameters are aj = l, A = l/5, and d = 1 (solid), d = 4 
(dashed), and d = 10 (dotted). 



decays to zero, while Pv{T) for the vertex model con- 
verges to a nonzero value for large d. 

The behavior of Pg (T) can be rather simply explained. 
Using the result of Eq. ([7]), for small A, 



P. 



,{T) « 2-"^ (l + « (1 - XT/4f « e-'^^^/4. 

. 

Thus, so long as A and T are independent of d, the hitting 
probability decays exponentially in d for the subspace 
model. 

The behavior of Pv{T) is somewhat more subtle. First, 
it is clear from Eq. (I38p that Ap„ < A, and thus e~^''"'* > 

that 



This can be used in Eqs. pd|) and (l42l) to show 



P.{T)^Y.9n{T)>e- 



AT 



(45) 



Tl = 



This lower limit is shown as the dotted curve in Fig. [31 
and is in good agreement with the limiting behavior of 
the actual Py{T). Most importantly, this limit on the 
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FIG. 3: The hitting probabilities Pv(T) (solid) and Ps(T) 
(dashed) as a function of hypercube dimension d. The pa- 
rameters are tj = 1, A = 1/5, and T — 7r/2. Also shown is the 
hmiting value exp(— AT) (dotted) as described in the text. 



hitting probability of the vertex model is independent of 
d. 

This striking distinction of the vertex and subspace 
models requires some discussion. First, as was noted after 
Eqs. ^ and PU|) . the subspace model appears to get 
less information (i.e. about subspaces not vertices) yet 
causes greater decay of the off-diagonal elements of the 
density matrix. This paradox can also be seen in terms 
of the Lindblad operators, of which the subspace model 
has 2d, while the vertex model has 2'^. Presumably more 
operators should cause more decoherence, but this is not 
the case. 

In fact, the subspace model involves more measure- 
ments than the vertex model. To understand why this 
is the case, note that both of these master equations can 
be derived from the following discrete-time process [8]: 

p{t + 6t) = {l~mp)Up{t)U^ +pJ2PjUp{t)U^Pj. (46) 

j 

This non-unitary evolution describes unitary evolution 
given hy U = e~^^^* interrupted by m measurements 
specified by the projectors Pj, and each measurement 
occurring with probability p — X5t. The results of these 
measurements are then ignored, leading to a loss of in- 
formation. The projectors here are equal to the Lindblad 
operators specified above. The operators for the subspace 
model can be grouped into m = d measurements (each 
with 2 projectors summing to the identity), while the 
vertex model has only one (m = 1) measurement (with 
2'^ projectors summing to the identity). Thus, while each 
measurement in the subspace model extracts less infor- 
mation, it actually performs d measurements, which in 
the end extracts the same amount of information as the 
vertex model {d bits). 

As a second and related point, the original work by 
Alagic and Russell used a decoherence rate X = p/d (note 
that in their notation d = n). Naively, this would remove 
the exponential decay seen above. However, they also 
used an energy scale to = k/d, such that the relevant pa- 
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FIG. 4: The entropy S{t) as a function of time for the quan- 
tum walk with decoherence in the vertex model. The pa- 
rameters are lu — 1, X — 1/5, and d = 1 (solid black), d — 4 
(dashed), and d = 10 (dotted), and each curve has been scaled 
by 1/d. For comparison the function 1 — e^"^' is also shown 
(solid gray). As described in the text, the entropy for the 
subspace model is identical to the d = 1 results (solid black). 



rameter AT — ttX/(2lu) = ■np/{2k) remains fixed, recov- 
ering the exponential decay of PsiT). Clearly, allowing 
the parameters of the model to vary with d can remove 
the exponential decay; the comparison presented above 
has the nice feature that in both the vertex and subspace 
models the system loses information at the same rate {d\ 
bits per second). 

As an alternative perspective, consider a physical hy- 
percube network of qubits designed to implement quan- 
tum state transfer f!^. In such a network, one expects 
that it is the number of operators in the system Hamil- 
tonian that depends on d, and not their strength. The 
difference between the subspace and vertex models then 
consists in how the qubits are coupled to an environ- 
ment (causing decoherence) . As shown in the Appendix, 
if each qubit has its own environment with an indepen- 
dent dephasing process, the corresponding master equa- 
tion coincides with the vertex model. However, it is also 
shown that if certain sets of qubits share multiple envi- 
ronments, chosen in such a way that qubits in the same 
subspace of the hypercube undergo a collective dephasing 
process, the corresponding master equation will describe 
the subspace model. Thus, the two decoherence models 
correspond to different physical systems. 

The heuristic information arguments given above are 
confirmed by numerical calculations of the entropy 
S{t) = — tr[p(i) log2 /o(i)], shown for the vertex model 
in Fig. m which has the qualitative form S{t) ^ d{l — 
e~^*) w dXt for small t. While these calculations are 
for the vertex model, in the exact solution of the sub- 
space model the total density matrix is equal to the ten- 
sor product of d copies of the d = 1 density matrix. Thus, 
the entropy for the d-dimensional case is simply d times 
that for d = 1, for which the vertex and subspace models 
are identical. It is also interesting to observe that the 
entropy production (dS/dt) is greatest when the quan- 
tum walk is halfway between the initial and final vertices 
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(e.g. t = T/2 = 7r/4). At this time the system would, 
in the absence of decoherence, be in a coherent superpo- 
sition over the entire hypercube. Such a state is greatly 
disturbed by the position (vertex or subspace) measure- 
ments, leading to a large increase in entropy. As seen 
in the figure, small oscillations in the entropy occur for 
longer times. These have twice the frequency of the os- 
cillations in the hitting probability, and occur each time 
the quantum walk is superposed over the hypercube. 

Finally, this comparison of information can be used 
to reconcile the vertex model with the exponential de- 
cay found by Kendon and Tregenna [llj . They observed 
that the discrete-time hitting probability Pd decayed as 
Pd ^ e^'^P, where p is the probability of making a mea- 
surement of the position (vertex) after each time step. 
Since the number of time steps for this quantum walk 
to hit the opposing vertex is ~ (i7r/2, the total amount 
of information (with d bits per measurement) extracted 
during the walk is approximately d^pTi/2. For a true 
comparison, then, I equate this information to the total 
information extracted in the two models above (~ d\T 
bits), and solve for p. The result, p « 2\T / (ird), leads to 
Pd ^ e^^^^^'^ , in approximate agreement with Eq. (US)) 
above. Thus, so long as one performs less than approxi- 
mately one complete position measurement for the entire 
walk, both the discrete and continuous vertex model will 
hit the opposite vertex with appreciable probability, in- 
dependent of d. 



V. CONCLUSION 

In conclusion, I have reconsidered the effects of de- 
coherence in the continuous-time quantum walk on the 
hypercube. By considering a particular model of deco- 
herence (the vertex model), I have found that a quantum 
state initially localized at one vertex will propagate to the 
opposite vertex with a probability bounded by a constant 
that is independent of the size of the hypercube. This is 
in agreement with an earlier study jllj of a discrete-time 
quantum walk, once the information lost to decoherence 
has been considered. This is in striking disagreement 
with another model of decoherence (the subspace model) 
[l3| . in which the same amount of information is lost, but 
in an apparently much noisier fashion. I now consider the 
interpretation of this result. 

The subspace model can be interpreted as extracting 
information by a sequence of d measurements, one for 
each dimension of the hypercube, whereas the vertex 
model uses a single measurement. By making more mea- 
surements to extract the same amount of information, the 
subspace model disturbs the quantum walk much more 
than the vertex model. The number of these measure- 
ments is proportional to the hypercube dimension, and 
each disturbance prevents the walk from reaching its des- 
tination. It is this proliferation of measurements that 
leads to an exponential decay of the hitting probability 
as the hypercube gets larger and larger. 



Further understanding this nontrivial relationship be- 
tween information, measurements, and decoherence in 
quantum walks may be of some value. Indeed, the role 
of decoherence in the mixing of quantum walks remains 
to be exploited in quantum algorithmsjSi] , and may have 
implications for biology and beyond 
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APPENDIX: QUANTUM STATE TRANSFER 

Quantum state transfer involves a network of qubits 
interacting with a Hamiltonian of the form 

H = i ^ hn^k (X.Xk + Y,Yk) (A.l) 
]<k 

where X, Y, and Z are the Pauli matrices for each qubit, 
and the coupling matrix is given by ^jk- Decoherence 
in quantum state transfer can be modeled by a master 
equation of the general Linblad form 

dtp = -t[n/h, p] + J2 >^AhpL] - \l]l,p - \pl]l,), 

.(A.2) 

with Lindblad operators Lj and rates A^. In this Ap- 
pendix I will illustrate how two different models of de- 
coherence for a hypercube network qubits lead to the 
quantum walk decoherence models discussed in the text. 

The first is an independent decoherence model, in 
which the Lindblad terms are given by 

E Up4 - ' +^ E (z^pz^ - p) ' 

(A.3) 

with = {Xj =F iYj)/2. This model involves indepen- 
dent energy decay (Ti, also called amplitude damping) 
and dephasing (T^, also called phase damping) processes 
for each qubit. Using this in Eq. (|A.2[) . and assuming 
that the initial density matrix p(t = 0) = |a;)(a:|, where 
\x) — Xx\Q • • • 0) {Xx is the Pauli operator for the qubit 
at location x), it will subsequently evolve to 

p(t)=po,o(t)|0)(0|+^p,,,(t)|a;)(y|, (A.4) 

x,y 

where |0) ~ |0 • • • 0). More general general initial states 
are considered in [lO|. These density matrix elements 
satisfy the differential equations: 

dtPo,o = Tf^ Y.X P^,-^ = ^r^(l - Po,o) 

^tPx,y " ^ '^^z {^x.zPz,y Px,z^z,y^ 

-T^^Sx.yPx.x - 2[(2ri)"i -I- r->,,y(i - ,5,,^). 

(A.5) 
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Note that Po,o + J2x Px,x = 1, and probability is decaying 
from the excited-state subspace to the ground state at 
rate l/Ti. To isolate decoherence within the excited- 
state subspace I let p^^y = e~^^'^^px,y, to find that px,y 
satisfies 

^tpx.y — i ^ ^^ {S^x,zPz,y ) "t~ Px.yi^ ^x,y\ 

(A.6) 

This will agree with the vertex model given by Eq.® if 
Hj k = ^j.k and A — ST^^. The effects of this type of 
decoherence on the fidelity of perfect state transfer are 
discussed in [lol |. 

The second model is a collective decoherence model, in 
which the Lindblad terms are given by 

^EE^^"' - \{Slo.,p}) , (A.7) 

with Lindblad operators 

{x\xj—a} 



with / the identity operator and the Pauli matrix 
for the qubit at location x. These operators have been 
chosen such that Sj^a\x) = —\x) if the j-th bit in the 
binary representation of x equals a, otherwise Sj^a\x) — 
+ \x). This choice is not unique: Sj^a = n{a;|a:j=a} 
is another valid choice. For both cases, using the same 
initial condition as before, and the result that 

%a|x)(y|%„ = {2Sx,^y, - l)\x){y\, (A.9) 

I find that po,o (t) — and the following differential equa- 
tion for px,y{t): 



This will agree with the subspace model given by Eq. pU]) 
if Hj^k — ^j.k and A = 2T7^. 
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